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Summary 

A capability has been developed that makes it 
possible to apply to measured flow field data the 
visualization tools developed to display numerical 
solutions for computational fluid dynamic problems. 
The measurement monitor surface (MMS) procedure was 
applied to the analysis of flow field measurements 
within a low-aspect-ratio transonic axial-flow fan 
rotor obtained with 2-dimensional laser anemometry. 
The procedure generates input for the visualization 
tools developed to display numerical solutions for 
computational fluid dynamic problems. The relative 
Mach number contour plots obtained by this method 
resemble the conventional contour plots obtained by 
more traditional methods. The results demonstrate 
that the MMS procedure can be used to generate 
input for the multidimensional processing and anal- 
ysis tools developed for data from numerical flow 
field simulations. They show that an experimenter 
can apply the MMS procedure to his data and then 
use an interactive graphics program to display 
scalar quantities like the Mach number by profiles, 
carpet plots, contour lines, and surfaces using 
various colors. In addition, flow directionality 
can be illustrated by display of vector fields and 
particle traces. The MMS procedure thus has the 
potential to produce better 

insight into turbomachinery and other types of 
fluid flow physics. 

Introduction 

Flow visualization is an important experimen- 
tal technique used to study complex physical phe- 
nomena. Data sets obtained are images, well suited 

Copyright © 1991 by the American Institute of Aeronautics 
and Astronautics, Inc. No copyright is asserted in the 
United States under Title 17, U.S. Code. The U.S. Govern- 
ment has a royalty-free license to exercise all rights under 
the copyright claimed herein for Governmental purposes. 

All other rights are reserved by the copyright owner. 


MMN. mmm 

COM EILOSTRATIQHS 



to 2-dimensional image processing and analysis. 

Flow visualization is also an important tool used 
to deal with the vast amount of data generated by 
numerical simulations. Numerical simulations 
produce multidimensional data that require addi- 
tional processing and analysis tools to display the 
information in a form that leads to a better under- 
standing and evaluation of simulated flows. A 
computational fluid dynamic researcher at a graph- 
ics workstation can display various portions of a 
computed flow field, observe the solution from 
different vantage points, and interact with the 
computed solution. 1-3 Having the ability to inter- 
act with the solution provides an opportunity to 
discover information that might not be anticipated. 

The vast amount of data generated by numeri- 
cal simulations matches the flow of data from 
experiments using advanced instrumentation such as 
laser anemometry. It would be advantageous if the 
multidimensional processing and analysis tools 
developed for data from numerical simulations were 
available to be used for multidimensional experi- 
mental data. Then an experimenter using an 
interactive graphics program could display scalar 
quantities like the Mach number by profiles, carpet 
plots, contour lines, and surfaces using various 
color intensities. In addition flow directionality 
could be illustrated by displaying vector fields 
and particle traces. This capability would provide 
the experimenter with the same opportunity as a . 
computational fluid dynamic flow researcher to 
discover information that might not be anticipated. 

This paper focuses on the application of a 
flow visualization scheme developed to aid the 
analysis of numerical data from fluid flow simula- 
tions to experimental data from laser anemometer 
measurements made in a transonic axial-flow fan 
rotor. 4,5 The measurements consisted of two compo- 
nents of velocity measured at many points in the 
fan flow velocity field. These measurements were 
made as part of an extensive experimental, theoret- 
ical, and computational turbomachinery research 
program being conducted at the NASA Lewis Research 
Center. 
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Several problems were encountered in the pro- 
cess of applying the flow visualization schemes 
used on results from fluid dynamic simulations to 
experimental data. One problem was the poor reso- 
lution of the fan blade location when the flow 
visualization was applied to the experimental data. 
A second problem was due to the relatively low 
density of the measurement points in the flow 
field. The problems were dealt with simultaneously 
by properly enriching the measured data set through 
the application of a measurement monitor surface 
(MMS) procedure. In the following sections , the 
MMS procedure will be described, and results of the 
procedure will be shown and discussed. 
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Nomenclature 

LA axial position, 2 

local speed of sound 

speed of sound based on plenum 
total temperature 

LA circumferential position, 6 

specific heat at constant (pressure 

distance between interpolation point and 
data point, Eg. (6) 

interpolated value at x,y 

measured value at point x fc 

Jacobian of transformation, Eq. (2) 

relative Mach number, Eq. (16) 

molecular weight of air, 0.02897 gm/(kgm 
mole) 

polytropic exponent 

grid cluster control functions, Eq. (1) 
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pressure 

flow field vector 

LA radial position, r 

universal gas constant, 8.314 
J/ (gm mole K) 

radial coordinate 

local temperature 

stagnation temperature 

plenum total temperature 

blade velocity, a function of r, Eq. (9) 

axial velocity 

tangential velocity 

absolute velocity magnitude, Eq. (14) 

magnitude of the relative velocity, 

Eq. (15) 

axial velocity component in blade 
frame of reference, Eq. (11) 

tangential velocity component in blade . 
frame of reference, Eq. (12) 

physical plane grid coordinates, Eq. (1) 

weighting function 

axial direction 

elliptic equation coefficients, Eq. (2) 

Equation (7) 

ratio of specific heats 
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change in tangential velocity, Eq. (10) 
tangential direction 

computational grid coordinates, Eq. (1) 

angular velocity in rpm 

angular velocity in rad/sec, Eq. (8) 


Subscripts 


index of measurement point 
local 
static 

axial component 
tangential component 
derivative with respect to 7J 
derivative with respect to £ 

Procedure 

The procedure discussed herein was developed 
for a specific experimental data set. The experi- 
mental data set chosen is reported in Refs. 4 
and 5. It consists of two components of velocity 
measured at many points in the flow field of a 
transonic axial-flow fan rotor. The measurements 
were made using a laser anemometer. While the 
procedure discussed herein can be used to generate 
input for any numerical flow field visualization 
methodology, the one used for development is PLOT3D 
which is described in Refs. 1 to 3. 

As a first step in developing the capability 
to apply to measured flow field data the visualiza- 
tion tools developed to display numerical flow 
field solutions, a simple flow physics model was 
used to generate density, two components of momen- 
tum, and energy at each point in the measurement 
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grid. The flow physics model was based on the 
premise that the flow field data were acquired 
using laser anemometer measurements along surfaces 
that approximate the location of the rotor design 
streamlines. The flow physics model is described 
in Appendix A. A description of the procedure for 
calculating the density, the two-dimensional momen- 
tum, and energy is given in Appendix B. 

In looking at the display of the results 
using data at the measurement points only, the 
first problem observed was that the measurements 
did not clearly define the axial-flow fan rotor 
geometry. The presence of a rotor was indicated by 
the absence of data. Data were taken at only 30 
positions in the axial projection. Only 10 posi- 
tions were inside the rotor. Fifty positions were 
used in the circumferential projection. Hence, the 
rotor was defined by the absence of data at only 
about 10 points in the axial direction, and thus 
the rotor appeared jagged and nonaerodynamic. 

The second problem was that the. density of 
measurement points in the axial direction was low 
so that the contour routine in the flow visualiza- 
tion program produced notched contours. Both of 
these problems were dealt with simultaneously by 
the measurement monitor surface (MMS) procedure 
which will be described next. 

In spite of these discouraging results, the 
simple flow physics model is the basis of the MMS 
procedure, and the first step in the MMS procedure 
is to use a simple flow. physics model to generate., 
the density, the 2-dimensional momentum, . and energy 
at each measurement point in a coprdinate frame 
fixed with the axial-flow fan rotor (Appendices A 
and B) . 

The procedure uses the following information 
from the experiment: 

1. Measurements of axial-flow fan rotor 

geometry „ . 
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2. Laser anemometry data that consists of the 
axial velocity, V^, and tangential velo- 
city, V Q , and the coordinates of all the 
measured points 

3. Test condition data (the axial-flow fan 
rotor velocity, inlet static pressure, 
inlet stagnation temperature) 

The second step in the procedure is to use a 
simple grid generation scheme to derive a monitor 
surface grid that includes the blade shape, employs 
streamline- type coordinates, and encompasses all 
the points in the measurement grid. The axial-flow 
fan rotor geometry and the measurement grid scheme 
information are used to generate the new grid. The 
grid generation scheme used will be discussed in 
the following section. The third step in the pro- 
cedure is the use of a scattered data interpolation 
scheme known as Shepard's method to interpolate the 
density, the two components of momentum, and the 
energy values from the measurement grid to the new 
monitor surface grid. The grid interpolation pro- 
cedure is discussed in more detail in the section 
on grid interpolation. The fourth step in the 
procedure is to use the visualization tools to 
display the monitor surface density, momentum and 
energy fields, and quantities that can be calcu- 
lated from them. In practice, the results are dis- 
played using PL0T3D on an Silicon Graphics IRIS 
workstation. 6 

Grid generation 

The requirements of the monitor surface grid 
were that it include the blade shape, employ 
streamline-type coordinates, and encompass all the 
points in the measurement grid. Many grid genera- 
tion techniques might meet these requirements. The 
grid generation procedure used in this investiga- 
tion is discussed in Refs. 7 to 10. 

The computer code used to generate the new 
grid employed the same logic as the code presented 
in Ref. 8. However, the FORTRAN code was adapted 
to run on a VAX-CLUSTER and modifications were made 
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in places like the data input and data output 
schemes and the size of the arrays. 


The grid generation scheme involves a mapping 
from a computational plane (£,77) to a physical 
plane (x,y). The mapping is constructed by speci- 
fying the desired grid points on the boundary of 
the physical domain. The grid is generated by 
solving a Poisson partial differential equation in 
the form 


Jf£ + J!£ 

dx 2 dy 2 

+ Za 

dx 2 dy 2 


1 ) 


(1) 


(where P and Q are known functions used to con- 
trol interior grid clustering. 

These equations are transformed to computa- 
tional space by interchanging the roles of the 
variable (£ f T)) and (x,y). 

This yields the following set of elliptic 
equations: 

ax u - 2jbr Ct) + yx^ = -j 2 (p Xi + Cx„) (; 

® - 2 Pnn + + oyj 

where 


s - + A 

p - + Wi 

y = A + yi 


j = 


d(.y, y) 


-Wi 




( 3 ) 


Numerical computations are performed in the 
uniform rectangular transformed plane, (£,77), to 
solve Eq. ( 2 ). The solution provides the (x,y) 
coordinates of each point in physical space. The 
axial-flow fan rotor geometry was used in a simply 
connected region mapping. For a simply connected 
region, Dirichlet boundary conditions can be used 
at all boundary points. 


8 



The measurement grid is 3-dimensional and is 
represented in radial, r (or RP), circumferential, 

6 (or CP), and axial, z (or AP), coordinates. The 
(r,0,z) coordinates are used in analytical expres- 
sions, while the (RP, CP, and AP) coordinates are 
used in describing the experimental results dis- 
cussed in Refs. 4 and 5. However, the data were 
principally collected in the cascade plane which is 
a plane defined by z (or AP) and 6 (or RP) coor- 
dinates for nearly constant r (or RP). The mea- 
surement grid thus obtained is unrolled by the 
transformation 


x = z 
y = r(z)6 


(4) 


The transformation (Eq. (4)) applies to the 
monitor surface grid so that the physical space 
(x,y) is actually (z,r (z)6) . In addition, it is 
used in the grid interpolation procedure that will 
be discussed next. - Consequently, here also the 
physical space (x,y) is actually (z,r(z )Q) . 

Use of Dirichlet boundary conditions means 
the solution must have prescribed values on the 
boundary. Considering a rectangular region in the 
( 77 , £) space, the boundary conditions are specified 
as follows. The inflow circumferential boundary 
condition is applied on the left side of the rec- 
tangle and the outflow circumferential boundary 
condition is applied on the right side. The blade 
suction surface is defined on the top of the rec- 
tangle and the blade pressure surface is defined on 
the bottom of the rectangle. The nonlinear simul- 
taneous algebraic equations (Eq. (2)) were solved 
iteratively using a successive over-relaxation 
technique (SOR). 

Grid Interpolation 

Interpolating the values of density, momen- 
tum, and energy from the measurement grid to the 
monitor surface grid is not trivial. Many schemes 
of curve and surface fitting are available. 11 
However, most of these schemes assume the data 
points lie on a complete rectangular lattice in the 
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xy— coordinate plane and function values (and possi- 
bly slopes) are available at each point. In the 
case under consideration, the points are not 
assumed to satisfy any particular conditions as to 
spacing and density. Consequently, they are 
referred to as being "scattered.* Scattered data 
interpolation schemes are reviewed in Refs. 12 
to 14. 

Because of its simplicity, the scattered data 
interpolation technique used was Shepard's meth- 
od. " This method is based on the premise that 
the value of a function F(x,y) at a point (x,y) in 
the monitor surface grid will be similar to the 
value of a point f^ in the measurement grid if 
the points (x,y) and (x k ,y k ) are close. Since the 
point in the monitor surface grid can be near many 
points in the measurement grid, the value of F(x,y) 
determined using Shepard's method is given by a 
weighted average as follows: 

If 

E ,y)f k 

F{x,y) = ^ (5) 

E w k( w >y) 

l C-l 


where the weights are a function of the distance 
between the points and decrease as the distance 
between points increases. As a result, measurement 
points at large distances from a monitor surface 
point do not influence its value significantly. 

In this Equation, F(x,y) represents the den- 
sity, momentum, or energy in the new coordinate 
system, f represents the density, momentum, or 
energy in the old coordinate system (measurement > 
grid) (x fc ,y k ), and N is the number of points used 
in the weighted average. 

The scattered data interpolation works as 
follows. Each point in the new grid is first 
located relative to a cell defined by the old grid. 
While it is preferable that the point in the new 
grid be in a cell defined by the old grid, it is 
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not necessary. This method has the distinct advan- 
tage that the point in the new gird need not be in 
a cell defined by the old grid. The four corners 
of this cell are used in the interpolation (N = 4). 
The procedure provided for a choice of weight, 

W k (x,y) , as either d k x ,d' k , or e' Bd *d* 2 where the 
distance between points is defined by 


d k - y/{ x ~ 


(r - Vkf 


and 



E 


(7) 


This value of a is selected to reduce the 
influence of far away points. - 

The first derivatives of F(x,y) vanish at the 
data points (x k ,y fc ). As a result, F(x,y) is flat 
near the data points. 

Shepard's method of scattered data interpola- 
tion was tried with each of the weights discussed 
in this section. No visual change was apparent in 
the results. 



To implement the measurement monitor surface 
(MMS) procedure one must write computer code. 
Consequently, one must have some experience with 
several disjoint engineering disciplines so that 
everything can be blended together. The MMS that 
makes it possible to apply to measured axial-flow 
fan rotor flow field data the visualization tools 
developed to display numerical solutions from com- 
putational fluid dynamic problems is based on a 
blend of knowledge from four fields. . The first 
component of the procedure is a simple flow physics 
model based on the premise that the flow field data 
were acquired using laser anemometer measurements 
along surfaces that approximate the location of the 
rotor design streamlines. Using this physical 
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model, information about the test condition, and 
the measured data, the procedure generates a den- 
sity, 2-dimensional momentum, and energy at each 
measurement point in a coordinate frame fixed with 
the axial-flow fan rotor. The second component of 
the procedure is a grid generation scheme that is 
used to derive the monitor surface grid that 
includes the blade shape, has streamline— type coor- 
dinates, and includes the points in the measurement 
grid (grid generation section). The third compo- 
nent of the procedure is the use of a scattered 
data interpolation scheme known as Shepard's method 
to interpolate the density, the two components of 
momentum, and the energy values from the measure- 
ment grid to the new monitor surface grid (grid 
interpolation section). The fourth component of 
the procedure is the visualization tools that are 
used to visualize the monitor surface and quanti- 
ties that can be calculated from them. 

In practice, the results were displayed using 
PLOT3D on an Silicon Graphics IRIS workstation. 6 
Consequently, some familiarity with the PLOT3D 
manual is also needed. In addition, one needs to 
understand certain aspects of the turbomachinery 
experiment such as the laser anemometry measurement 
scheme and the axial-flow fan geometry. 

The flow field calculations and the use of 
PLOT3D are discussed in Appendices A and B. The 
flow logic of the computer programs developed to 
implement the procedure is discussed in Appendix C. 

Results 


In this section, the measurement monitor 
surface procedure will be demonstrated. In addi- 
tion, a relative Mach number (Appendix A) plot 
obtained using the MMS procedure and the PLOT3D 
program will be compared to a plot of the same data 
made using conventional methods. The experiment 
will be discussed first. Then results obtained 
with the MMS procedure and the PLOT3D flow visual- 
ization program will be presented. 
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The Experiment 


Although the PL0T3D program makes it possible 
to view many variables, in this paper the emphasis 
is on the shock structure. This is displayed in 
plots of the relative Mach number (Appendix A) . 

The shock structure to be shown is from measure- 
ments acquired in a low-aspect-ratio transonic 
axial-flow fan rotor. Much of the total pressure 
rise which occurs in advanced fans and core com- 
pressor inlet stages occurs across the rotor pas- 
sage shock. Consequently, these measurements were 
made to help develop accurate models of the passage 
shock that could be used in the axial-flow fan 
rotor design process. Details of the experiment 
and the results of the test program are discussed 
in Refs. 4 and 5. The rotor aspect ratio is 1.56 
and the design tip relative Mach number is 1.38. 

The flow field shown is for a near-stall operating 
condition. The plots are based on intrablade ve- 
locity measurements acquired with a laser fringe 
anemometer on blade-to-blade planes in the super- 
sonic region at 30-percent span. 

A schematic diagram of the NASA Lewis single- 
stage axial-flow compressor test facility is shown 
in Fig. 1. The facility is discussed in Ref. 4. 

The test rotor, NASA rotor 67, is shown in Fig. 2. 
The rotor is also discussed in Ref. 4. 

The laser anemometer (LA) measurement 
location in the stationary reference frame is spec- 
ified by the parameters AP, RP, and CP. As shown 
schematically in Fig. 3, AP is the axial position, 
RP is the radial position measured from the com- 
pressor axis of rotation, and CP the circumferen- 
tial position measured relative to the z-0-plane. 
Velocity measurements occurred along the circum- 
ferential measurement line which is swept through 
the LA measurement location as the blades rotate 
past the measurement position. The 2-dimensional 
surface of interest is the (x,y) plane which corre- 
sponds to the (z,r(z ) 0 ) plane. Detailed informa- 
tion on the LA measurements is given in Ref. 4. 
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Laser anemometer axial measurement locations 
(AP) and radial measurement locations (RP) are 
shown in the meridional view of the test fan rotor 
presented in Fig. 4. The laser anemometry results 
presented were obtained at the 30-percent span 
location shown in Fig. 4. Radial surveys of total 
and static pressure, total temperature, and flow 
angle were made at survey stations 1 and 2 located 
upstream and downstream of the rotor as shown in 
Fig. 4. 


Figure 5 from Ref. 4 shows a contour plot of 
the relative Mach number at a near-stall flow test 
condition at 30-percent span from the shroud pro- 
duced by a nonworkstation graphical procedure. The 
MMS procedure was applied to the data measured at 
this same test condition. Since the same relative 
Mach number calculation procedure is used and the 
same data were used, it expected that the MMS pro- 
cedure results will be similar to the ones shown in 
Fig. 5. 

The MMS Procedure 


Figure 6 shows the grid used for the laser 
anemometer velocity measurements. The grid has 
30 horizontal divisions and 100 vertical divisions. 
Note that the axial-fiow fan rotor shape defined by 
the absence of data in the measurement grid has a 
jagged appearance. The grid generated by the MMS 
procedure for this example is 51 by 102 and it is 
shown in Fig. 7. 

Figure 8 shows the measured relative Mach 
number in the measurement grid. Each shade of 
color represents a region of constant Mach number. 
Note the jagged appearance of the axial-flow fan 
rotor which is an artifact of the measurement 
scheme and a significant distracting influence. 

In addition, the low density of points (30) in the 
axial direction may be hindering the contour draw- 
ing package. After viewing these results, the MMS 
procedure was proposed as a possible solution to 
the problems shown in this figure. 
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The MMS procedure results are shown in 
Fig. 9. Note that the axial-flow fan rotor shape 
is now smooth and the contours are smoother. Fig- 
ures 10 and 11 show closeups of the measured and 
MMS procedure-derived relative Mach number con- 
tours. As expected, the MMS procedure results are 
similar to those presented in Fig. 5. 

Looking at the plot of the grid produced by 
the MMS procedure used herein which was made using 
the PLOT3D flow visualization program, the shape of 
the ends of the axial-flow fan rotor are not shown 
as being as sharp as one expects. The problem is 
related to the boundary conditions specified in the 
(£'*?) Pl ane * The (x,y) coordinates along the blade 
pressure surface are specified on the bottom £ 
axis as boundary conditions. The (x,y) coordinates 
along the blade suction surface are specified on 
the top £ axis as boundary conditions. The tips 
of the blade were specified by three points on the 
pressure surface and three points on the suction 
surface. The points specifying the ends are close 
together while those along the blade are spaced 
much further apart. This spacing problem caused a 
great amount of grid distortion at the blade ends. 
In addition, the distortion was enhanced by the 
PLOT3D visualization program. While the grid used 
was adequate for demonstrating the MMS procedure, 
it is obvious that further refinement is necessary. 
A more detailed investigation of the grid distor- 
tion problem at the ends of the blade would lead to. 
resolution of this problem. However, it might 
require using a second grid program to adjust the 
grid points at the blade tips. 

A different type of distortion occurs in the 
flow field at the downstream ends of the axial-flow 
fan rotor. Here the elongated low Mach number 
region at the end of the rotor in Fig. 8 is now a 
sequence of oval low Mach regions in Fig. 9. The 
shape of the elongated low Mach number region is 
attributed to a lack of data in this region. The 
sequence of oval low Mach regions is an artifact of 
the scheme which tries to interpolate the measure- 
ments onto a fine grid. 
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While only relative Mach results are shown, 
all the standard 2-dimensional flow quantity calcu- 
lations that the PLOT 3D program can make are avail- 
able. These results thus indicate that an experi- 
menter can apply the MMS procedure to his data and 
then use an interactive graphics program to display 
scalar quantities like the Mach number by profiles, 
carpet plots, contour lines, and surfaces using 
various colors. In addition, flow directionality 
can be illustrated by display of vector fields and 
particle traces. The MMS procedure thus has the 
potential to produce better insight into turbo- 
machinery and other types of flow physics. 

Discussion 


Graphical techniques are used by experiment- 
ers and computational fluid dynamic researchers to 
provide an understanding of flow field phenomena 
and the effects of changes in configuration. In 
addition, computational fluid dynamic researchers 
use graphic displays to assess the accuracy of a 
calculation, and to evaluate errors, inconsisten- 
cies, or improper modeling. 

The results presented herein suggest experi- 
menters could interactively find a new use for 
graphical methods since they could assess the 
usefulness of a particular measurement grid in 
obtaining data for code validation by using the 
techniques presented herein to move computational 
flow field predictions into a measurement grid. 
Viewing computational results can indicate regions 
of large gradients where more points are desirable. 
Simulated experimental results for a particular 
measurement grid scheme would supply additional 
information on the the shape of the contours for 
physical quantities such as relative Mach number. 

This technique could also be used during an 
experimental research program to allow experiment- 
ers to monitor interactively their results in a 
display scheme that included the experimental geom- 
etry while the data were being collected. 
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Computational fluid dynamic researchers often 
need to use overlapping grids. The results pre- 
sented herein suggest they could use Shepard's 
method to interpolate between these grids. This 
might be especially useful for grids adjusted by 
direct interaction or grids adjusted during compu- 
tation. 


Conclusion 


A capability has been devloped that makes it 
possible to apply to measured flow field data the 
visualization tools developed to display numerical 
solutions from computational fluid dynamic 
problems. The procedure is based on a blend of 
material from four fields of knowledge. First, a 
simple flow physics model is used to generate den- 
sity, two components of momentum, and energy at 
each point in the measurement grid. The flow phys- 
ics model is based on the premise that the flow 
field data were acquired using laser anemometer 
measurements along surfaces which approximate the 
location of the rotor design streamlines. Second, 
a simple grid generation scheme is used to derive a 
grid which includes the axial-flow fan rotor shape 
and streamline-type coordinates for all the points 
in the measurement grid. Third, a scattered data 
interpolation scheme known as Shepard's method is 
used to move information from the measurement grid 
to the new grid forming a monitor surface for den- 
sity, the two components of momentum, and energy. 
Finally, tools developed to visualize numerical 
fluid flow results are used to visualize the moni- 
tor surfaces and quantities that can be calculated 
from them. 

The measurement monitor surface procedure 
(MMS procedure) has been demonstrated using laser 
anemometry measurements in a transonic axial-flow 
fan rotor made at a near-stall flow test condition 
at a position 30-percent of the span from the 
shroud. The Mach number results obtained by the 
MMS procedure for this case are similar to results 
obtained using standard techniques. 
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These results indicate that a MMS procedure 
can be used to generate input for the multi- 
dimensional processing and analysis tools developed 
for data from numerical flow field simulations. 

They show that an experimenter can apply an MMS 
procedure to his data and then use an interactive 
graphics program to display scalar quantities like 
the Mach number by profiles, carpet plots, contour 
lines, and surfaces using various colors. In addi- 
tion, flow directionality can be illustrated by 
display of vector fields and particle trances. The 
MMS procedure thus has the potential to produce 
better insight into turbomachinery and other types 
of flow physics. 

Appendix A - Flow Calculations 

A procedure for putting laser anemometer (LA) 
flow field measurements in a transonic axial-flow 
fan rotor into a form suitable for display in the 
measurement grid by PLOT 3D is discussed. 

One limitation of the procedure is that it is 
assumed that the results from the LA measurements 
have been acquired along measurement surfaces which 
approximate the location of the rotor design 
streamlines. The approximation is clearly incor- 
rect when the grid points are inside a blade sur- 
face. 


In the analysis of rotor flow and energy 
transfer, the variables of concern. are of interest 
in a coordinate system rotating with the blade. In 
this section the relationships used to convert the 
measured axial and tangential velocity to the ro- 
tating coordinate system will be set forth. 

Basic thermodynamic equations used herein 
such as definitions of specific heat, the speed of 
sound, the perfect gas law, work done by the sys- 
tem, and the energy equaton are given in thermody- 
namic text books such as Ref. 15 and in discussions 
of turbine design (Refs. 16 and 17) and these 
expressions are not repeated. 
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The velocity vector relationships used are 
presented here since the data are shown in a frame 
of reference rotating with the blade while it is 
measured in a stationary laboratory frame of refer- 
ence. However, turbine design textbooks such as 
Refs. 16 or 17 should be consulted for more 
information. 

Velocity Vectors 

The angular velocity in rad/sec is related to 
the velocity in rpm as follows: 


| co | = 2rcQ/60 (8) 

where 

| Q | is the velocity in rpm. 

For this case |Q| = 16042 rpm 
and |b/| = 1679.9 rad/sec. 

For the system under consideration the fol- 
lowing conventions are used. The angular velocity 
is taken to be negative. 

The blade velocity is 

U(r) =-r(z) | <o | . (9) 

For this case with r = 21 cm, the blade 
velocity is -352 m/s. 

The change in tangential velocity is 

AV 0 = (V e ) ini9t -(V e ) wit (10) 

The laser measures two components of an abso- 
lute velocity. An axial component, .V , and a tan- 
gential component V^. For flow in and across the 
stators, the absolute velocites are of interst . 

For flow in and across the rotors, the veloc- 
ities are considered relative to the velocity of 


19 



the, rotating blades. In terms of relative veloci- 
ties, flow in a rotating blade row can be analyzed 
in a manner similar to the analysis of flow in a 
stationary passage. 

Since the blade velocity is always in the 
tangential direction, the axial velocities are 
unchanged. 

W z = V z ( 11 ) 

In the coordinate system used the tangential 
velocity at the rotor entrance (-100 percent Chord) 
is zero. At the rotor exit (-157.23 percent Chord) 
the tangential velocity is -155 m/s. The relative 
tangential velocity is 

W 9 = V, + r |o»| ( 12 ) 


Relative Mach Number 


Using basic thermodynamic equations for the 
specific heat, the speed of sound, the work done by 
the system,' and the energy equation, a quantity 
called the local speed of sound, a^ , is defined as 
follows: 



(Y - 1) 


^-Av^rco 



(13) 


where 


= (Vi ♦ V?) ,/2 ( 14 ) 

and a Q is the speed of sound based on the plenum 
total temperature, r is the radius, 7 is the 
ratio of the specific heat at constant pressure to 
the specific heat at constant volume, W is the 
rotational speed, and Av^ is the change in tan- 
gential velocity from the plenum to the point of 
interest. This local speed is only an approxima- 
tion since radial velocity components are neglected 
in calculating the absolute velocity. 
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The magnitude of the relative velocity is 

W = (wl + W%) 1/2 ' ( 15 ) 


The relative Mach number is the given by 

M r .1 = -£ ( 16 ) 

a l 

Appendix B - Using PLOT3D 


Measurement grid 


The measurement grid is 3-dimensional and is 
represented in radial, r, circumferential, 6, and 
axial, z, coordinates. However, the data was prin- 
cipally collected in the cascade plane which is a 
plane defined by z and 6 coordinates for nearly 
constant r. The measurement grid thus obtained is 
unrolled by the transformation 


y = r(z)Q 


(17) 


Flow Field 

To use PLOT3D to analyze the results from a 
computational fluid dynamic simulation one must 
have as input to the PLOT3D program the computation 
grid and values at each grid point of the flow 
field vector Q where the components of Q are 
the density (Q x = p) t the three compoents of momen- 
tum (Q 2 - pV x rQ 3 = /0V ,Q4 = pV z ) , and the stagnation 
energy per unit volume (Q 5 = pc^ T q ) . 

Since PLOT3D is to be used to visualize 
2-dimensional information in a frame of reference 
rotating with the blade, the flow field vector, ~Q, 
must be specified in terms of the relative veloci- 
ties . 


The mapping used to specify Q was chosen 
based on Eq. (17) which was used to define the 
grid. 
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The following relations are used to specify 
these quantities in the rotating blade coordinate 
system: 


a. = P, 

Q 2 " P W z 

a, - P« t 


a, = o.o 



+ 


W 7. 

2 


where the relative velocity W, is given by 

Eq. (15), the relative tangential velocity, W 0 , is 

given by Eq. (12), the local temperature, , is 


T « - 


2 

M w a < 

~w 


N 


the local acoustic velocity, a^ , is given by 
Eq. (13), 7 is the ratio of specific heats, R is 
the universal gas constant, 8.314 J/(gm mole K), 
and M # is the molecular weight of air, 0.02897 
gm/ (kgm mole) . 

Using basic thermodynamic equations for the 
specific heat, the speed of sound, the work done by 
the system, and the energy equation, we find the 
local stagnation temperature to be given by 

c A T *) t s c e T t + = c p T 0 - AV 0 r« (18) 


or 


( T *) t = 


= T. + 


II IZl 
2 \ Y 


\ 5* 
I R 


' abs 



R 


AV 9 r<* 

(19) 
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The stagnation temperature, (T s )^ , is also 
called the total temperature. It is the tempera-" 
ture attained when a gas at static local tempera- 
ture, , and velocity, V abg , are brought to rest 
adiabat Ically . 


The measured initial static pressure, 

(P s ) i n i e t' and the stagnation temperature, (T 8 ) ±nlet , 


defined above 

are used as 

initial conditions. 

We ' 

assume that 

inlet 

" <Ps> C 

(20) 


inlet 

= < T s>C 

(21) 


The measured initial pressure, (p ) , is 

inlet 

regarded at the pressure of a fluid brought to rest 
isentropically from a velocity V abg and a static 
local pressure p^ . Thus 


Pt 

(P B ) Inlet 


( T s) 


inlet 


■ T 
(T-l) 


( 22 ) 


From static temperature, T^ , and static pres- 
sure, p^ , the density p ^ is calculated using the 

ideal gas equation. 

Between points along a streamline the flow 
process is considered to be a polytropic process so 
that the pressure and temperature at these points 
are related as follows: 

££ 1 (23) 

Again at each local point, using T^ and 
p^ , the density p. is calculated with the ideal 
gas equation. This local density is used to deter- 
mine the input for PLOT3D. 

In addition to the calculations of the PLOT3D 
input, the total or stagnation pressure is 
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calculated at each local point assuming an isentro- 
pic process: 


(P«) 


( T s\ 




■ 1 
(r-i) 


(24) 


The pressure, temperature, and density calcu- 
lated from the laser measurements are an approxima- 
tion to the true pressure, temperature, and density 
since only two of the three components of velocity 
are available. Thus, the Mach number is also an 
approximation to the true Mach number. 

Appendix C - Software 

Computer programs were developed or were 
available to perform each of the four functions. 

The prime impediment in implementing the MMS proce- 
dure was the need to combine these programs in a 
single system. For the benefit of anyone who must 
implement a MMS procedure the system logic which 
was developed will be discussed. In this discus- 
sion the names of the data files and Fortran 
programs will be used. Figure C-l shows a logic 
diagram of the files and three FORTRAN software 
programs used to put the procedure into practice. 
Note that two of the three programs are related to 
grid generation. The input data files contain the 
axial-flow fan rotor geometry information 
(SURF 10 . DAT ) and the 2-dimensional laser anemometry 
information (ST_NS30.DAT and SD_NS30.DAT). 

Three software programs are used. The first 
program ( INFLAP01 . FOR) computes, using spline curve 
fits if necessary, the input boundary information 
for the grid generation program (FOIL2.FOR). It 
produces files with status information that are 
used to study the behavior of the program 
( INFLAP01 .OUT) , the grid information needed by the 
third program (GRID1.DAT) , and input information 
for the second program which generates the monitor 
surface grid (AIRF0IL5.DAT) . 


24 



The second program (FOIL2.FOR) generates the 
monitor surface grid (GRID1 . X, Gridl .Q) using a 
successive over-relaxation (SOR) to solve the par- 
tial differential equation that, defines the grid. 

In addition, it produces files with status informa- 
tion that indicate the rate of convergence of the 
SOR iteration procedure (AIRFOILS .OUT) and the 
initial grid guess (INOUT1.X, INOUTl.Q). . The moni- 
tor surface data sets GRID1.Q and INOUTl.Q are 
empty. However, the (IN0UT1.X, INOUTl.Q) and 
(GRID1.X, GRID1.Q) file pairs can be viewed using 
PLOT3D to see if the new grid is satisfactory and 
assess the starting grid. 

The third program (COMBINE2 . FOR) combines the 
measured grid information, the measured 
2-dimensional velocities, the physical model, the 
experiment test conditions, and the new grid, to 
produce the monitor surfaces of density 
2-dimensional momentum, and energy (BINDEA.X, 
BINDEA.Q) using Shepard's method to interpolate 
between the monitor surface grid and the measure- 
ment grid. This program also produces monitor 
surfaces on the original measured grid ( BINDED.X, 
BINDED.Q) and status information ( COMBINE. OUT ) . 

Each of these programs has a specific data 
file with control information (INFLAP01.DAT, 
FOIL2.DAT, and COMBINE2.DAT). In addition, there 
is a master data file with control information 
(MASTER.DAT). These control files perform func- 
tions such as specifying the number of points in 
the monitor surface grid, the grid cluster control 
parameters, and the weighting procedure to be used 
by Shepard's method. 
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FIG. C-1 . - BLOCK DIAGRAM OF COMPUTER PROGRAM LOGIC. 
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FIG. 1. - NASA LEWIS SINGLE-STAGE COMPRESSOR TEST FACILITY. 



ORIFICE 


INLET THROTTLE VALVES-*— 


DRIVE 

MOTOR 


EXHAUST 


CD-10916-11 


STRAIGHTENER 
TEST SCREEN - 
STAGE - 


GEAR- 


PLENUM - 


/ 

/ 

/ 

/ 

'-COLLECTOR/ 

THROTTLE 

VALVE 


ALTITUDE 

EXHAUST 


ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 



FIG. 2. - TEST ROTOR. 


LA OPTICAL AXISn 

i 



FIG. 3. - DEFINITION OF LASER ANEMOMETER CIRCUM- 
FERENTIAL MEASUREMENT LINE. 
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FIG. 4. - MERIDIONAL VIEW OF TEST FAN ROTOR SHOWING 
LASER ANEMOMETER AND AERODYNAMIC SURVEY LOCATIONS. 



FIG. 5. - CONTOUR PLOT OF RELATIVE MACH NUMBER AT NEAR- 
STALL FLOW TEST CONDITION AT 30-PERCENT OF THE SPAN 
FROM THE SHROUD. 
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FIG. 6. - MEASUREMENT GRID. 



FIG. 7. - MONITOR SURFACE GRID. 
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CONTOUR t.EUELS 



FIG. 8. - SHADED SURFACE RELATIVE MACH 
NUMBER PLOT BASED ON MEASURED DATA. 



FIG. 9. - SHADED SURFACE RELATIVE MACH 
NUMBER PLOT BASED ON MONITOR SURFACES. 
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FIG. 10. - CLOSEUP VIEW OF SHADED SURFACE RELATIVE MACH NUMBER PLOT BASED ON MEASURED DATA. 
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FIG. 11. - CLOSEUP VIEW OF SHADED SURFACE RELATIVE MACH NUMBER PLOT BASED ON MONITOR SURFACES. 





